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Coulomb-gauge ghost and gluon propagators in SU{3) lattice Yang-Mills theory 
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We study the momentum dependence of the ghost propagator and of the space and time compo- 
nents of the gluon propagator at equal time in pure 5*17(3) lattice Coulomb gauge theory carrying out 
a joint analysis of data collected independently at RCNP Osaka and Humboldt University Berlin. 
We focus on the scaling behavior of these propagators at /3 = 5.8, . . . , 6.2 and apply a matching 
technique to relate the data for the different lattice cutoffs. Thereby, lattice artifacts are found to 
be rather strong for both instantaneous gluon propagators at large momentum. As a byproduct we 
obtain the respective lattice scale dependences a(/3) for the transversal gluon and the ghost propa- 
gator which indeed run faster with /3 than two-loop running, but slightly slower than what is known 
from the Necco-Sommer analysis of the heavy quark potential. The abnormal a(/3) dependence as 
determined from the instantaneous time-time gluon propagator, -D44, remains a problem, though. 
The role of residual gauge-fixing influencing Dn is discussed. 
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I. INTRODUCTION 



Lattice investigations of the gluon and ghost propa- 
gator have become an important topic over the last ten 
years after the pioneering lattice studies in the Landau 
gauge appeared in the late eighties and nineties [1, 2, 3, 
4, 5, 6, 7] and after the coupled solutions to the corre- 
sponding Dyson-Schwinger equations in the deep infrared 
momentum region were found [8, 9]. Since then, the 
available amount of lattice data on these propagators has 
grown (see, e.g., [10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20]) 
and also studies based on functional methods have made 
considerable progress [21, 22, 23, 24] such that in Lan- 
dau gauge one is nowadays in the comfortable situation 
to confront continuum results with a broad set of inde- 
pendent lattice data (see, e.g., [25, 26] for recent discus- 
sions). 

In Coulomb gauge, comparably few lattice investiga- 
tions of the aforementioned propagators have been per- 
formed. For example, Langfeld and Moyaerts [27] as 
well as Cucchieri and Zwanziger [28], and very recently 
also Burgio, Quandt and Reinhardt [29, 30] have carried 
out such computations for the gauge group SU{2). In 
fact, the Coulomb gauge provides an interesting alter- 
native to the Landau gauge since the resulting Hamilto- 
nian approach allows to apply the variational principle 
to get analytic results for the QCD vacuum wave func- 
tional and for the spectrum of hadronic bound states. 
This approach has been mainly pursued by the Tiibingen 
group [31] in recent years. Their investigations of (trun- 



cated) systems of Dyson-Schwinger equations in the 
Coulomb gauge provided various solutions in the in- 
frared, including also asymptotic "conformal" or "scal- 
ing" solutions that are characterized in the infrared by 
a singular power-law behavior of the ghost propagator 
and a powerlike vanishing transversal gluon propaga- 
tor [32, 33, 34], very similar to what was found in Landau 
gauge. 

Among us, the authors from Japan have done sev- 
eral lattice investigations before, concentrating on the 
Coulomb gauge for 5[/(3) gauge fields. The instanta- 
neous gluon propagators and the ghost propagator were 
computed in [35] , and correlators of incomplete Polyakov- 
loops were also determined. The latter was studied in or- 
der to interpolate between the confinement potential Vc 
(derived from Wilson loops) and the Coulomb potential 
Vcoui (known to restrict Vc from above [36, 37]). Further- 
more, it was possible in this way to extract potentials for 
the quark-antiquark singlet and octet channels, as well as 
for the quark-quark symmetric sextet and antisymmetric 
anti-triplet channels [38, 39]. The eigenvalue spectrum 
of the Coulomb-gauge Faddeev-Popov (FP) operator was 
studied in [40, 41]. 

Recently, some of us have also computed the gluon and 
ghost propagators as well as the Coulomb potential Vcoui, 
the latter directly from the FP operator, however [42]. 
For Vcoui very strong Gribov-copy effects were reported, 
and it still remains difhcult to give a final answer for the 
infrared momentum limit [43]. Independent of that, the 
factorization assumption proposed in [44] relating Vcoui 
to the square of the ghost propagator was found to be 
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strongly violated at low momenta. 

In this paper we present a joint analysis of data from 
Berlin and Osaka for the instantaneous propagators of 
both ghost and gluons. For the transverse gluon propa- 
gator as well as for the ghost propagator similar infrared 
properties as for the Landau gauge are expected also 
for the Coulomb gauge. In Landau gauge, e.g., a gluon 
propagator vanishing in the zero-momentum limit (or a 
infrared-diverging ghost dressing function) is crucial from 
the point of view of the Gribov-Zwanziger confinement 
scenario [45, 46] or the Kugo-Ojima confinement crite- 
rion [47]. In Coulomb gauge, the instantaneous time-time 
gluon correlator should become singular and be related 
to the effective Coulomb potential. 

This study provides a comprehensive set of lattice data 
on the instantaneous gluon and ghost propagators in the 
Coulomb gauge of pure SU{3) lattice gauge theory. We 
discuss their momentum dependence and analyze in de- 
tail apparent scaling violations of the space and time 
components of the gluon propagator. Wc show that these 
violations can be ameliorated if different cuts are applied 
on the data. In fact, they are effectively eliminated by a 
matching procedure that provides us also with the run- 
ning of the lattice scale a{(3), separately for each propa- 
gator. With the exception of the time-time propagator 
£'44, we find these runnings to be in good agreement with 
other prescriptions. The behavior of the renormalization 
coefficients, that are also provided by the matching pro- 
cedure, is smooth a long as /3 > 6.0. 

The structure of the paper is as follows: In Sect. II 
we describe the setup of our lattice simulation including 
details on our gauge-fixing algorithms. Sect. Ill intro- 
duces the relevant lattice observables. The data for the 
propagators is discussed in Sect. IV where we report on 
obvious scaling violations for the gluon propagators. We 
then use a matching procedure to relate the propagators 
for different lattice cutoffs to each other and discuss the 
outcome of this for the instantaneous gluon propagators 
and the ghost propagator in Sects. V, VI and VII. We 
present in detail the interplay of the matching procedure 
with the necessity of an additional momentum cutoff that 
restricts the reliability of the data to relatively small mo- 
menta \pa\ < a. The lattice scale dependence a(/?), as 
determined thereby, is compared to what is known from 
the literature in Sect. VIII. Finally, in Sect. IX, we dis- 
cuss the momentum dependence of the propagators in 
both the ultraviolet and infrared region. Wc draw our 
conclusions in Sect. X. To make the paper self-consistent, 
we give a brief outline of the matching procedure for 
the Coulomb-gauge propagators in Appendix A. Fit and 
matching tables are presented in Appendix B. 



II. LATTICE FIELD ENSEMBLES AND GAUGE 
FIXING 

The results discussed below are based on an extensive 
set of quenched gauge configurations generated in Os- 



aka and Berlin. At both places we employed Wilson's 
one-plaquette action and a standard heatbath algorithm 
(including microcanonical steps) for thermalization, but 
used different values of the inverse coupling /? and dif- 
ferent lattice sizes L^. Those, together with a couple of 
other useful parameters, are listed in Table I that can be 
found in Appendix B. 

In our analysis below we combine the data from Osaka 
with the data obtained in Berlin. Both sets are nicely 
consistent with each other as we checked by comparing 
data at (3 = 5.8 and 6.0. 

Configurations were fixed to Coulomb gauge via max- 
imizations of the Coulomb gauge functional 



^£/[5] = EE^^^TrC/f(f,i), 



(1) 



i=l S.t 



where Uf{x,t) = g{x,t) Ui{x,t) g'' {x + i,t). For a fixed 
U this was done by iteratively changing g using a stan- 
dard overrelaxation (OR) algorithm in Osaka and the 
simulated annealing method combined with subsequent 
overrelaxation (SA-fOR) as in Refs. [42, 43] in Berlin. 
Strictly speaking, different gauge-fixing methods may 
cause variations in the data of gauge- variant observables 
due to the Gribov ambiguity. This is, in particular, true 
for the ghost propagator at very low momentum (with 
deviations up to 5%, see Ref. [43] for a detailed account 
on that). 

Each maximum of Fu [g] automatically satisfies the lat- 
tice Coulomb gauge condition 







(2) 



for all color components (c — 1,...,8). Here is 
the lattice backward derivative in one of the three spa- 
tial directions i, and Af is the lattice gluon field. Via 
Af, = A^^T", with the eight generators T" of SU{3) in 
the fundamental representation, the gauge field compo- 
nents arc defined in terms of the gauge-fixed links C/^ {x, t) 
through 



A„ 



2iago 



U^{x, t)-Ul{x, t) 



tracclcss 



(3) 



where go is the bare coupling (related to (3 — Q/ g^) and a 
denotes the lattice spacing. Note that we follow the mid- 
point definition which defines A^ at the midpoint of a link 
U^{x, t), i.e., Ai = Ai{x + ^i, t) and A4 = A4{x, t + i). 

Obviously, maximization of i^;7[,9] proceeds indepen- 
dently in each time slice, as neither the Coulomb gauge 
functional (1) nor the resulting gauge condition (2) fixes 
a link in temporal direction. We observe that the time 
slices of a given configuration may behave very differently 
during the iterative gauge-fixing process. In fact, we find 
that the number of necessary iterations may differ by a 
factor of 10 to 20 between the individual time slices of 
a given configuration. This obstruction reficcts that a 
topological tunneling might happen within one or a few 
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subsequent time slices. In some cases there were time 
slices which could not be fixed within a certain prede- 
fined number of iterations. Then, the gauge-fixing pro- 
cess was repeated for those time slices starting from a dif- 
ferent randomly chosen gauge transformation g restricted 
to that slice, while leaving the "well-behaved" (already 
gauge-fixed) time slices untouched. In the majority of 
cases, time slices did not show any recalcitrancy during 
gauge fixing, though. 

After all the individual time slices were maximized, the 
original configuration U was gauge-transformed 

Ui{x,t) g{x,t)Ui{x,t)g\x + i,t) , (4a) 
U4ix,t) g{x,t)U4ix,t)g'<{x,t + l) , (4b) 

i.e., including also the time- like links. 

After having fixed the Coulomb gauge by maximizing 
the functional ^[/[s] there is still freedom to carry out 
gauge transformations ht G SU (3) which only depend on 
time. One way to fix this residual gauge freedom is to 
maximize the functional 



^f/4(f,t) 



(5) 

where U is the Coulomb-gauge transformed configura- 
tion (Eqs. 4). Links in time direction arc finally gauge- 
transformed under h(t) as 

Ui{x,t) ^ h(t)Ui{x,t)h^t + l) (6a) 
whereas spatial links are transformed as 

U^{x,t) ^ h{t)U,{x,t)h^t) , (6b) 

which preserves the Coulomb gauge. Equal-time observ- 
ablcs involving only spatial links, i.e., the transversal 
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FIG. 2: Mean link values {^9\eTrU^) for difTerent direc- 
tions n at 13 = 6.0 as function of the linear lattice size L. 
Notice the strong volume dependence of (|$HeTrf74) after 
residual gauge-fixing (RGF). 



gluon or the ghost propagator, are not affected by this 
residual gauge freedom. 

If the transversal or the time-time gluon propagator 
was to be defined for non-equal time, the residual gauge 
would have to be fixed as well. For equal times, however, 
it is not clear to what extent, if at all, the results for the 
instantaneous time-time gluon propagator DfKif) would 
change if the remaining gauge freedom was fixed. We 
will check this in Sect. VI by comparing data for for 
fixed residual gauge freedom (Berlin data) to that where 
this freedom was left unfixed (Osaka data). 

In Fig. 1 we show the average trace (|9leTr[/j) for 
spatial links after Coulomb gauge fixing (invariant under 
residual gauge fixing) and also the average (|5KeTr?74) 
before and after residual gauge fixing as function of /3. 
The data is taken for a 18^ lattice. Without residual 
gauge-fixing the average trace of time-like links vanishes, 
whereas after residual gauge-fixing the expectation value 
is finite and increases with increasing /3. Fig. 2 demon- 
strates that the average (^ Die Tr U4) after residual gauge- 
fixing (shown for f3 = 6.0) steeply decreases with increas- 
ing lattice volume, in contrast to (i 9^e Tr Uj) for spatial 
links. 

In Sect. VI we will demonstrate that the difference 
of D44 measured with and without residual gauge-fixing 
can be completely accounted for by a multiplicative, 
momentum- independent rescaling, e.g., by normalizing 
the matched D44 propagator from both versions at some 
reference scale p = ^J.. Apart from this, the residual gauge 
fixing has an impact only on the value of the propagator 
at zero momentum but this is not of importance for our 
present study. 



FIG. 1: Mean link values (ifReTrC/^) for different direc- 
tions /i on a 18"* lattice. Within errors the spatial links 
($ReTrf/i>, i = 2,3 (not shown) are equal to {JHeTrC/i>. 
(| SHe Tr U4) refers to the time- like links before and after resid- 
ual gauge-fixing (RGF). 



III. COULOMB-GAUGE PROPAGATORS ON 
THE LATTICE 

The space and time components of the gluon field eval- 
uated in momentum space enter the bare instantaneous 



4 



transversal and time-time gluon propagator as the Monte 
Carlo correlators 

D'l^'ik) = (i^(fc)i5(-fc)) , (7a) 

Dttik) = (i^(fc)ii(-fc)) . (7b) 

Here Ai and denote the spatial Fourier transforms 
of the lattice gluon fields at a fixed time t with integer 
momenta ki € {—L/2, L/2]. An average over all time 
slices is understood. Dfj' is diagonal in color space and 
transverse in momentum space. On the lattice it takes 
the form 

Dt^k) ^ S'^" - P^^^MM.^ D,,^p) (8) 

with 

K(fc)^^sin(^^) . (9) 

This is simply due to the lattice Coulomb gauge condition 
which in momentum space translates into 

3 

Y,p,{k)A,ik) - (10) 

i=l 

for all k. In the following we use p = | p{k) | to simplify 
the notation wherever applicable. 



When analyzing data on D°'j it is natural to associate 
the physical momentum with p. Lattice results then re- 
produce the continuum tensor structure of D°^. Devia- 
tions from its tree-level form are described by the dimen- 
sionless dressing function Ztr(p), defined by 

Ztr=pDtr{p). (11) 

Analogously, the time-time gluon propagator Df| may 
be presented in the form of either D44(p) or ^44(73). Both 
are related to the full propagator through 

D2t{k) = S''" Diiip) - £!lM . (12) 

P 

Depending on the particular focus, data below is pre- 
sented in either one or the other form. 

The ghost propagator is defined as the expectation 
value of the inverse Faddeev-Popov (FP) operator M 

{{M-^X^^)^5u.5''''G{x-x') (13) 

at a fixed time t = t' (subsequently averaged over all 
time-slices). The FP operator is local in time and, by 
virtue of the chosen Coulomb gauge functional (1), has 
on the lattice the dimensionlcss form 



Mf^, = Sw J2 Tr [ {T% T''} (u,ix, t) + U,{x - 1 1)) 6^3. 



2 T'T^ U,ix, t) 5,^., 2 T'^T<'U.{x ~ z, t) 5,_:^., 
I 



(14) 



where i is a unit vector in spatial direction, x = {x, t) 
and T° is a generator of SU (3) in the fundamental rep- 
resentation. 

We are particularly interested in the momentum de- 
pendence of the ghost dressing function 



J{p') = (p{k)) G{p{k)), 



where 



'M- 



(15) 



(16) 



Working in momentum space, it is convenient to invert 
M for a selection of momenta and colors c forming right- 
hand side plane-wave sources £,'^{k) = go-Cg^^rik-x/L 



k = {k,0) ^ 0. We use a preconditioned conjugate- 
gradient algorithm described in Ref. [13], adapted to 
Coulomb gauge, to accelerate the inversion of M . Alter- 
natively, we could have used a selection of point sources 
^ and Fourier-transformed the vectors \M^^^{x) provid- 
ing an estimator for the ghost propagator at once for 
all momenta, however with less statistical accuracy. The 
plane- wave method automatically ensures that J{j?) is 
averaged over all time slices. Moreover, translational in- 
variance is exploited to improve the estimator. Note that 
M cannot be inverted for fc = due to its eight trivial 
(constant) zero eigenmodes. 

Multiplicative renormalizability is a well-established 
property of the gluon and ghost propagators in Lan- 
dau gauge. In Coulomb gauge, to our knowledge, this 
has been proven yet only up to one-loop by Watson and 
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Reinhardt quite recently [48]. Their result for the bare 
dressing functions obtained in the 4D momentum space, 
formally translated from dimensional to lattice regular- 
ization looks as follows (omitting possible lattice correc- 
tions) 

Z^^ = 1 + 5o Ctr [log {a\x + y)) + H^O] + 0{gt), 
= l+gl C44 [log {a'ix + v)) + i/44(C)] + 0(5^), 
= 1 + gl C.J [\og{a^y) + H.j] + 0{gt), (17) 

where the momentum variables are x = pi,y — p'^,C — 
x/y, and the lattice cutoff is a~^. Ctr, C44, Cj and Hj 
denote constants. Note the non-trivial dependence on 
C for both the Coulomb-gauge gluon dressing functions. 
When multiplicatively rcnormalizing the dressing func- 
tions, e.g., in a momentum subtraction (MOM) scheme 
at some scale this dependence has to be carefully taken 
into account (sec also [30]), in particular for equal-time 
correlators which according to Eq. (17) require an inte- 
gration over p4 or (. Of course, 3> /i has to be 
ensured for that which, admittedly, is very difhcult to 
achieve in nowadays lattice computations. Even more, 
lattice computations are typically carried out at several 
values of (3, i.e., at different cutoff- values. In general one 
should expect, that the corresponding dressing functions 
at different lattice spacings, say a and a, arc related to 
each other by a finite renormalization of the Z-factors 
which will obviously depend only on the ratio a/a and 
not on the momenta. This will then hold also for the 
dressing functions Eq. (17) and correspondingly also for 
the equal- time correlators. 

In the case of Landau gauge those Z-factors turned out 
to be close to unity for the gluon and ghost dressing func- 
tions at similar values of /?. Therefore, it is more or less 
sufficient for them to express the various lattice spacings 
by a unique physical scale, e.g., via the Sommer-scalc 
parameter ro = 0.5 fm and the interpolation formula of 
Necco and Sommer [49] 

ln(a/ro) = -1.6804 - 1.7331(/3 - 6.0) 

+ 0.7849(/? - 6.0)2 _ o.4428(/3 - 6.0)^ (18) 

obtained from the lattice analysis of the static quark- 
antiquark potential and applicable in the range 5.7 < 
f3 < 6.92. The remaining lattice artifacts were sufficiently 
dealt with by applying cone and cylinder cuts to the mo- 
menta [50]. While a cone cut addresses finite- volume ef- 
fects, the cylinder cut is an easy and effective method to 
reduce artifacts due to the broken rotational symmetry. 
We shall apply both these cuts also to our data shown 
below. However, in what follows we will demonstrate 
that the approach, even if sufficient for Landau gauge, 
is not quite enough for the case of Coulomb gauge. In 
fact, beside applying the usual cone and cylinder cuts 
one has to restrict momentum components to also satisfy 
api < a < 2 and to apply non-trivial finite renormaliza- 
tions between the different cutoff values. 

Also, we shall not use a priori the Necco-Sommcr scal- 
ing relation [Eq. (18)] but instead find the specific scaling 



behavior for each of the propagators defined above and 
present their data in terms of the finest available lattice 
scale at either (3 = 6.20 or 6.10. For this we employ the 
matching procedure of Ref. [50] adapted here to Coulomb 
gauge. A detailed outline of this method applied to our 
propagators is given in Appendix A. 



IV. DISCRETIZATION ERRORS OF THE BARE 
LATTICE DATA 

We start our discussion by revisiting the strong scaling 
violations we reported for the transversal and the time- 
time gluon propagator in [35, 42]. There, we used the 
interpolation formula Eq. (18) to assign physical units to 
the lattice momenta and applied a multiplicative normal- 
ization at ^ = 2 GeV for all values of p. This procedure, 
however, leads to serious disretization errors for both the 
instantaneous transversal and the time-time gluon prop- 
agator (see Fig. 3), whereas the ghost propagator looks 
much more satisfactory in this respect (see Fig. 4). 

Challenged by these scaling violations, in the next sec- 
tions we shall perform a matching procedure that merges 
the data for different f3 into one bare lattice propagator 
associated with the highest available lattice cutoff. 

In a first step, however, we consider here the scaling 
violations and argue them to indicate that the admissi- 
ble range of lattice momenta needs to be restricted even 
further than what the cylinder and cone cuts would do. 
For this, we introduce a new momentum cut that will 
be applied in addition to those two cuts. Basically, not 
the full Brillouin zone should be eligible when analyz- 
ing the corresponding propagator data, but only that at 
momenta (Eq. (9)) whose components are restricted to 
\pia\ < a < 2. For the sake of brevity we will refer to 
this cut as the "a-cut" in what follows. 

Fig. 5 illustrates the effect of the a-cut on the instan- 
taneous transversal gluon propagator. Note that in this 
figure (as in [35, 42]) we have used the Necco-Sommer 
formula (18) (and tq = 0.5 fm) to assign physical units 
to momenta and propagator. Obviously, when decreas- 
ing a less and less data points survive this cut but those 
that do show a much better overlap than before (see, in 
particular, the lower panels of Fig. 5). 



V. MATCHING THE TRANSVERSAL GLUON 
PROPAGATOR 

Still, the disagreement between data from different P 
does not completely disappear. Therefore, in a next step, 
we relax the a priori universal a{f3) dependence (e.g., 
that according to Eq. (18)) and apply the matching pro- 
cedure of Ref. [50] as explained in Appendix A. It pro- 
vides us with multiplicative renormalization factors de- 
pending on the ratios of the lattice spacings and with 
the specific dependence of the lattice spacing a — a((3) 
separately for each propagator. 
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FIG. 3: The instantaneous transverse gluon propagator (left) and the dressing function of the instantaneous time-time gluon 
propagator (right) using the Necco-Sommer scaling relation and normalized sd, jj, — 2 GeV. The data refers to approximately 
equal physical volumes and has been cylinder and cone cut. Data was produced at HU Berlin. 
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FIG. 4: The ghost dressing function using the Necco-Sommer 
scaling relation and normalized at jj, — 2 GeV. The data has 
been cylinder and cone cut. Data at /3 = 5.8, 6.0, 6.2 was 
collected at HU Berlin, data at /3 = 5.9, 6.1 at RCNP Osaka. 



We start with the instantaneous transversal gluon 
propagator and first match data obtained on two lattices 
with approximately the same physical volume, i.e., data 
on a L** = 32'*, /? = 5.8 lattice with data on a = 48^, 
(3 ~ 6.0 lattice. Besides the cone and cylinder cuts we 
apply two different a-cuts (with a = 0.6 and a = 0.5) 
before performing the matching procedure. Our aim is 
to compare the influence of the a-cut on the quality of 
matching. The result, with L"* = 32"*, (3 = 5.8 being the 
coarse and = 48^, (3 = 6.0 being the fine lattice, is 
shown in Fig. 6 for both a-cuts. 

We obtain good matching of both data sets with a bet- 
ter result for a = 0.5 (see the /dof listed in Table II). 
There is hardly any difference between the best result 
of the matching procedure on one hand and directly im- 
posing the Necco-Sommer scaling relation on the other 
(Fig. 5). Indeed, our matching procedure nearly repro- 
duces the lattice-spacing ratios as given through Eq. (18) 



(see Table II) in Appendix B. 

Next we extend the matching to all values of /3 = 
5.8 . . . 6.2 using data obtained on 32^ and 48^ lattices. 
Since (3 = 6.2 has the finest lattice spacing, the matching 
is performed between data at /? = 6.2 (setting the refer- 
ence scale) and data at all other /3 = 5.8 . . . 6.1. We also 
compare the result for four different a-cuts. The results 
are summarized in Table III. 

We not only find the ratios of lattice spacings to rise 
monotonously upon decreasing /3, but also the ratios 
of the renormalization constants to be about the same 
(somewhat below 1.0), i.e., almost independent on (3. 

As a general rule, smaller values for a result in lower 

/dof values and hence better matching. We also see 
that the matching procedure nearly always results in 
a lattice-spacing ratio smaller than that given through 
Eq. (18), although the discrepance decreases with a taken 
smaller. As shown in Fig. 7 for the best a-cut (a = 0.5), 
we achieve a virtually perfect matching of the instan- 
taneous transversal gluon propagator over all data ob- 
tained at /3 = 5.8 . . .6.2. Comparing with the /dof for 
a = 1.0, we conclude that applying the a-cut is essential 
to achieve a good overlap of the data. Our combined fi- 
nal result indicates a flattening of the propagator in the 
infrared region which is worth to be explored further. We 
expect a tendency to show a plateau as recently seen in 
the Landau gauge case [51], which excludes a vanishing 
gluon propagator in the infrared limit. 



VI. MATCHING THE TIME- TIME GLUON 
PROPAGATOR 

We now apply the matching procedure to the instan- 
taneous time-time gluon propagator 1)44. As for the 
transversal propagator Dt^.^ we first match data obtained 
on two lattices with approximately the same physical vol- 
ume, i.e., data for = 32*, (3 = 5.8 with data for 
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FIG. 5: Effect of the a-cut on the bare lattice data of the instantaneous transverse gluon propagator. The data (all for a 32* 
lattice) is presented imposing different «-cuts: no a-cut (top left), \pia\ < 1 (top right), \pia\ < 0.6 (bottom left), \pia\ < 0.5 
(bottom right). Cylinder and cone cuts have been applied as before. Necco-Sommer scaling (with ro = 0.5 fm) has been used 
to get a(/?) and physical momenta p. The data for 13 = 5.8, 6.0, 6.2 (5.9, 6.1) was obtained at HU Berlin (RCNP Osaka). 
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FIG. 7; The instantaneous transverse gluon propagator ob- 
tained by matching all data from 32^ and 48* lattices, in- 
cluding data from both Osaka and Berlin collected at five 
different P values. The result is shown for a fixed a-cut with 
\pia\ < 0.5. Data are cylinder and cone cut. 

L"* ~ 48** and (3 = 6.0, respectively. We also compare 
the effect of two difFcrcnt a-cuts (a = 0.6 and a — 0.5) 
in addition to the usual cylinder and cone cuts. 

Whereas the matching seems to work reliably as shown 
in Fig. 8 and as demonstrated in Table IV by the low 
X^/dof value for a = 0.5, the obtained lattice spacing 
ratio is now significantly larger than predicted by the 
Necco-Sommer scaling relation (see Table IV). This is in 
striking contrast to what we have observed in the case of 
the transversal gluon propagator. 

We now merge all data for D44 in the interval /3 = 
5.8 . . . 6.2 obtained on 32^ and 48^^ lattices. However, 
since the instantaneous time-time gluon propagator is 
more sensitive to Gribov copies [42] and since we have 
employed different gauge fixing procedures at HU Berlin 
and at RNCP Osaka, we first match the corresponding 
data sets separately. The resulting fit parameters are 
summarized in Table V (Osaka) and Table VI (Berlin). 

Matching the Osaka data one finds that the ratios Ra 
of lattice spacings rise monotonously upon decreasing /3 
but much stronger than in Eq. (18). The ratio of renor- 
malization constants is still compatible with unity for 
(3 = 6.0 if compared to /3 = 6.1 (providing the reference 
scale), but it decreases abruptly between /3 = 6.0 and 
(3 = 5.9. The x^/dof value is acceptable only for an 
Qf-cut where a = 0.5. 

The Berlin data allows only to compare (3 ~ 6.0 and 
(3 = 5.8 to /3 = 6.2 (which sets the reference scale). The 
ratios of the lattice spacings are compatible with the re- 
sults for the Osaka data. The ratio of the renormaliza- 
tion constants is still compatible with unity for (3 = 6.0, 
if compared to (3 = 6.2, but drops between /3 = 6.0 and 
(3 = 5.8 similar to the Osaka data. The x^ /dof is unac- 
ceptably large. 

In Fig. 9 we show our final result for the instantaneous 
time-time gluon propagator having matched and com- 
bined all the Osaka and Berlin data. We have now used 



a unique scale set by a(/3 = 6.0) to give all momenta in 
physical units. The quality of the fits in this case is worse 
compared with the fits of the transverse gluon propaga- 
tor. Nevertheless, the scaling behavior does look quite 
reasonable and there is an improvement compared to the 
results presented in Table IV. The reason is probably 
that we have moved closer to the continuum limit by in- 
cluding data from [3 = 6.1 and (3 = 6.2. In the infrared 
region the data points obtained in Berlin and Osaka split. 
We interprete this as a consequence of the use of differ- 
ent gauge fixing techniques. The more efficient simulated 
annealing method weakens the singular behavior as seen 
also for the ghost propagator [43]. 

In passing, we revisit the question whether there is 
a difference between the instantaneous time-time gluon 
propagator if one is applying the residual gauge-fixing 
(Berlin data) or not (Osaka data). Fig. 10 shows the re- 
sults for (3 = 6.0. The propagator at momentum pa ^ 
seems not to depend on the volume, but on the proce- 
dure (cf. the left panel). The latter is understandable if 
one looks back at Fig. 1 and there at the difference for 
the time-like links. It is remarkable that the difference 
between the two cases can be eliminated by a uniform 
multiplicative rescaling. This is accomplished by nor- 
malizing the propagator I?44/a at pa = 2.0 to 1.0 (cf. the 
right panel). The residual gauge fixing has only an im- 
pact on the value of the propagator at zero momentum, 
D44{ap = 0)/a. With residual gauge fixing this value is 
obviously smaller as expected. 



VII. MATCHING THE GHOST PROPAGATOR 

Finally we apply the matching procedure to the ghost 
propagator. As we have seen in Fig. 4, Necco-Sommer 
scaling is only weakly violated. Therefore, we expect a 
matching result that closely follows this behavior. 

The fitting results (respectively for the Osaka and 
Berlin data) are presented in Tabs. VII and VIII, there 
relative to the highest /3max = 6.1 and /3max = 6.2 in both 
cases. For no a-cut (a = 2.0) the ratios of lattice spacing 
a(/?)/a(/3inax) reproduce almost perfectly the correspond- 
ing ratios according to Eq. (18). Nevertheless, we notice 
that a smaller a leads to some deviation from that scal- 
ing, in particular at the lowest /3. Including the results 
of the separate fits, the Osaka and Berlin data are after- 
wards combined in Fig. 11 showing there the result only 
for the most restrictive a-cut (a = 0.5). 

Matching the Osaka data yields an overall very good 
x'^/dof. The ratios of the renormalization constants are 
all compatible with unity, and the ratios of lattice spac- 
ings rise monotonously upon decreasing (3. When no 
a-cut is applied the fitted Ra rise in accordance with 
Eq. (18), while restricted a-cuts lead to R^s which grow 
slightly slower. Matching the Berlin data results in the 
same tendencies, but the x^/dof turned out to be very 
large. Probably due to the fact that the Berlin ghost 
data are averaged over all time slices and the Osaka data 
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FIG. 8: The instantaneous time-time gluon dressing function after matching Berhn data from two lattices of approximately 
equal physical volume and applying cylinder and cone cuts. Two a-cuts are compared: \pia\ < 0.6 (left) and \pia\ < 0.5 (right). 
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FIG. 9: The dressing function of the instantaneous time-time 
gluon propagator after matching Berlin and Osaka data sep- 
arately on 32'* and 48* lattices for five different /3 values. Fit 
parameters underlying the matching are listed in Tabs. V and 
VI. Data has been cylinder and cone cut and results are shown 
for a fixed a-cut with \pia\ < 0.5. The two fits for the mo- 
mentum dependence in the IR are described in Sect. IX. 



are not, the errors of the Berlin data are smaller by an 
order of magnitude. 



VIII. THE SCALING BEHAVIOR OF THE 
PROPAGATORS 

At this point wc can check now if our individual re- 
sults on a(/3) reproduce a unique running lattice scale. 
We had started with Necco-Sommcr scaling, but aban- 
doned this, fully relying on the matching procedure to 
produce the "correct" lattice scale function. Let us re- 
mind that the lattice scales as found may deviate from 
asymptotic scaling (which of course is strictly valid only 
for (3 — > 00) and also from that derived for other observ- 
ables (e.g., Necco-Sommcr scaling derived for the static 



quark-antiquark force in pure SU{3) gauge theory). 

With Fig. 12 we summarize our results on the fitted 
scaling behavior of the lattice spacing in terms of the 
finest lattice spacing available in our study. For all the 
propagators considered we have plotted the ratios of the 
lattice spacings relative to the finest lattice (at /? = 6.2 
for the transverse gluon propagator in the left panel, and 
at /3 = 6.1 for 1)44 and the ghost propagators in the 
right panel, respectively) as a function of /3 as found 
through the matching procedure with two choices of the 
a-cut. For the transversal gluon propagator the data 
points (corresponding to both choices of a-cuts) fit very 
well between the curves corresponding to two-loop run- 
ning of the lattice spacing, 

and to the relation (Eq. (18)). The same is observed for 
the ghost propagator. In this case, also applying no cut 
(a = 2.0) results in a reasonable result. 

On the other hand, the a(/?) dependence as found 
from the matching procedure of the instantaneous time- 
time gluon propagator is not only much stronger than in 
Eq. (18), but also inconsistent with the scaling law for the 
other propagators. We cannot say to which extent such 
a faster running is beyond some general bound and un- 
fortunately have to conclude that the problem of the bad 
scaling behavior for the instantaneous time-time gluon 
propagator remains unsolved yet. 



IX. FITTING THE BEHAVIOR AT LARGE AND 
SMALL MOMENTUM 

Having successfully merged data for the propagators 
from simulations at different /3 values, one may try to fit 
their ultraviolet (UV) behavior and partly also to extract 
some infrared (IR) exponents. 
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FIG. 10: Left: comparison between the unrenormalized instantaneous time-time gluon propagator D44 with and without the 
residual gauge fixing. Right: the same but formaUy normahzed at pa = 2. The cyhnder cut has apphed. Berhn and Osaka 
data for /3 = 6.0 and lattice sizes 32'' and 48'' are shown together. 
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FIG. 11: The ghost dressing function after matching data 
from 32'' and 48* lattices including data from both Osaka and 
Berlin at five different /3 values. For the fits see Table VII and 
Table VIII. The result is shown for a fixed a-cut: \pia\ < 0.5. 
Both the cylinder and the cone cut are applied as usual. 



For the transverse gluon propagator we try a power- 
law ansatz 



1 



Vtr 



(20) 



to describe the behavior at large momenta. For the UV 
fitting we use the data points above some minimal mo- 
mentum (|p|min in units oia~^{P = 6.2)) and investigate 
the dependence of the anomalous dimension r]tr on the fit- 
ting range and the a-cut. Fit results are collected in Ta- 
ble IX and our best fits (with a = 0.5) give Tytr = 0.40(2). 
Qualitatively, the behavior we find is similar to the UV fit 
given in Ref. [27], though there [for SU{2)] a somewhat 
bigger exponent r^tr ~ 0.5(1) was found. 

For the longitudinal gluon dressing function Z44 we try 



(21) 

The results are collected in Tabs. X and XI, respectively. 
This dressing function was not studied in Ref. [27]. 

For the ghost dressing function, analogous to Ref. [27], 
we adopt a logarithmic ansatz in the UV region 



mr = 



Cgh 



ln(|j^/Acoui)^ ' 
and a power-law ansatz for the IR behavior 



(22) 



(23) 



Fits results for either momentum region are given in 
Tabs. XII and XIII, respectively. The UV fits scatter 
with a (the additional momentum cut), though, the most 
stable results are obtained for no a-cut [a — 2.0). With 
a suitably restricted fit interval fits are stable and give 
Acouia(/3 = 6.0) = 0.275(20) or Acoui^o 1.37(10) and 
7 = 0.33(1). For SU{2) this exponent was found to be 
7 = 0.26(2) [27]. 

The IR fits are quite stable and give Kgh = 0.435(6) 
without applying a-cuts, even though we admit that the 

/dof values are rather large. In Ref. [27] a value Kgh = 
0.49(1) was found (corresponding to 2k there). 



X. CONCLUSIONS 

We have investigated the momentum dependence of 
the instantaneous ghost and gluon propagators of pure 
SU{3) lattice Coulomb gauge theory. Our study repre- 
sents a joint analysis of data from lattice simulations in- 
dependently performed at Berlin and Osaka for the Wil- 
son gauge action in the range /3 = 5.8, . . . , 6.2. 
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FIG. 12: The scaling relation as it emerges from matching the data obtained in simulations with the Wilson action. Left: for 
the transversal gluon propagator with two a-cuts for comparison. Right: for the ghost propagator with no a-cut {a = 2.0) and 
an a-cut with a = 0.6 for comparison, and for the time-time gluon propagator with an a-cut with a — 0.6, a rather conservative 
momentum cut. 



For these values of /3, we find apparent scaling viola- 
tions for both the spatially transversal and the time-time 
gluon propagator, while for the ghost propagator such 
violations are surprisingly mild. Our inspection of the 
gluon propagator data shows that the violations there 
are basically due to data that survives a cylinder cut but 
involves momentum components close to the upper end 
of the Brillouin zone. Consequently, if additionally an 
a-cut like \pia\ < 0.5 is applied to the data, scaling vio- 
lations are under much better control. The price to pay 
are strong restrictions of allowed momenta which, in our 
opinion, should not only satisfy the cylinder and cone 
cuts but also \pia\ < 0.5 (a-cut). This is the first result 
of our paper. 

Second, we find that the scaling violations can be suf- 
ficiently reduced if, in addition to the aformentioned 
cuts, a matching procedure (see Appendix A) is used 
to merge data. That is, instead of imposing one partic- 
ular a{(3) dependence (e.g., that of Rcf. [49]) and nor- 
malizing the data for the different lattice cutoffs such 
that they coincide at a particular reference scale, both 
the a(/3) dependence and the relative normalization fac- 
tors are determined through an optimization method 
that seeks the best overlap of data. It turns out, that 
the matching procedure applied to either the transversal 
gluon or the ghost propagator provides us with a a(/3) 
dependence only slightly different from what is known 
from [49], somewhere in between Necco-Sommer scaling 
and asymptotic two-loop scaling. Note that the matching 
procedure would allow us to fix the lattice spacing if we 
were to simulate also beyond the interval 5.7 < (3 < 6.92 
covered by the Necco-Sommer analysis. 

Generally we can say that the matching analysis re- 
sults in ratios of the renormalization constants closer to 
unity at P > 6.0. Future lattice studies of gluon and 



ghost propagators should be performed in that region. 
The fact that - except for the ghost propagator - the 
matching performs better the more restrictive a-cuts are 
applied shows that the momenta with components close 
to the upper end of the Brillouin zone are far from the 
continuum limit. This might signal a more general effect, 
namely that observables closer to the infrared region have 
better scaling properties. 

Unfortunately, we could not correct the scaling viola- 
tions for the instantaneous time-time gluon propagator. 
For this, these violations are so strong that the a(/3) de- 
pendence as found through the matching is far from what 
we find for the other propagators. In fact, a(/3) in this 
case is found running too fast. Moreover, for (3 < 5.9 the 
ratio of renormalization constants drops compared to the 
behavior at (3 > 6.0 such that the assumptions and re- 
sults of the matching analysis for the propagator 
must be considered with caution. 

We mention that for the SU (2) transversal gluon prop- 
agator it has been argued [30] that the correct instan- 
taneous propagator can be reconstructed only from the 
full 4-dimensional space-time propagator. There, a resid- 
ual gauge-fixing was applied that enforces A4 = const. 
Therefore, it needs to be scrutinized whether the scal- 
ing violations, that we have seen here for the transver- 
sal gluon propagator, are really due to the alleged 
(multiplicative) non-renormalizability of the Coulomb 
gauge [52] when residual gauge-fixing is applied or not. 
Our results for the transversal propagator suggest a more 
mundane resolution: exclude too large momenta from the 
analysis and allow for an independently determined run- 
ning lattice spacing, then data within a very restricted 
range of momenta (in units a) can be successfully merged 
and gives a a(/3) dependence that agrees with what is 
known form the literature. 
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We stress again that our result for the time-time gluon 
propagator is non-acceptable. The a(/3) dependence as 
found for this is far from the running scale for the other 
propagators. In the light of this, the argument of non- 
renormalizability might still be valid for the A4 compo- 
nent of the gluon field. 

When applying fits to the data at cither low or large 
momenta (though restricted by quite stringent bounds) 
we obtain qualitatively similar UV and IR fits as reported 
for the SU{2) theory in [27]. 
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APPENDIX A: MATCHING PROCEDURE 

In this appendix we describe the matching procedure 
of [50] applied to Coulomb gauge. The procedure does 
not rely on any given lattice scale dependence a(/3) but 
allows us to extract this for each propagator individually. 

Under the assumption that the fixed-time gluon and 
ghost propagators in Coulomb gauge can be renormalizcd 



multiplicativcly (see Sect. HI and Ref. [48]), we aim at 
an optimal overlap of bare propagator data from a coarse 
lattice (with unknown lattice spacing Oc) and a fine lat- 
tice (with a lattice spacing aj that might be known). 
Using the fact that the bare and dimcnsionless lattice 
propagator Z?^ is a function of the product of the three- 
momentum p with the lattice spacing a only (the depen- 
dence on P is of course kept in mind) , and assuming that 
multiplicative renormalization is valid, the bare propaga- 
tors on the fine and coarse lattice arc related by 

afDf{paf) = Rz{af/ac) ■ acD^{pac) (Al) 

The renormalization factor Rz only depends on the ratio 
of the lattice cutoffs Ra ~ af/uc- Taking the logarithm 
gives 

\nDj;{paf) = \nD^{pac) - InRa+lnRz . (A2) 

Expressing the momentum on the coarse lattice in terms 
of the momentum on the fine lattice by 

flc = In(pac) = ln(pa/) - In i?a , (A3) 

we arrive at 

\nDf[ln{paf)] = In Df- [In (pa/) - In i?a] - Ini?^ + \nRz 
= In IHpaf) + Aa] + Az , (A4) 

where Ra = e"^" and Rz = 6-^"+'^^. 

Notice that and Az arc positive. We find the values 
for Ra and Rz from a fitting procedure as follows. 

Suppose that we have one data set {x = paf, Dj, <Jj}i 
with i ^ 1, . . . ,nf for the fine lattice and one data set 
{y = pUc, , <Jc]j with j = 1, . . . , rtc for the coarse lat- 
tice with a denoting the statistical error of the propaga- 
tor , and n/ resp. ric denoting the number of data 
points for the propagator on the fine and coarse lattice, 
respectively. Then, we use a fit to optimally match 
both data sets, i.e., to find the optimal overlap of the 
bare lattice propagator from the fine and the coarse lat- 
tice. To be specific, we minimize 



E 



D^x,) 



Rz 

Ra 



D. 



L int 



(ft) 



E 



Ra 

Rz 



jjL int 



(A5) 



In the first term is represented by the measured val- 
ues at the momenta pi (expressed as function of Xi — 
Pittf) and the corresponding error (Tf_,;, while is 
evaluated at these momenta by a cubic spline interpola- 
tion of the data for . In the second term, the role 



of Df" and Dcj^ is interchanged with respect to genuine 
data {uj = pjQc in -Dct) and interpolation of Df". With 
this definition of the matching is done as follows: 

1. Vary Aq over an interval (0, 1] with step size 0.001 
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FIG. 13: Details of matching Z*"^ measured on a L* = 32", (3 = 5.8 lattice with the data obtained on a L" = 48*, /3 = 6.0 
lattice. An a-cut with \pia\ < 0.5 was applied before matching. Left: ^ jdof as a function of the ratio of lattice spacings 
Ra. = ac(5.8)/a/(6.0). Right: the 68.3% confidence region spanned by ratios of lattice spacing Ra — ac(5.8)/a/(6.0) and 
renormalization constants Rz = .^(5.8)/Z(6.0) determined through the matching procedure. 



and determine the optimal Az giving the lowest 
/dof for each value of considered, 

2. Identify the best overall combination of Aa and Az 
by searching for the global minimum of / dof . 

This provides us with the optimal choice of Ra and Rz- 
The error of Ra and Rz are given by the 68.3% confidence 
region, i.e., the region of fit parameters Ra and Rz with 

/dof < Xmin/^o/ + 1. An illustration of this is given in 
Fig. 13 for matching the instantaneous transversal gluon 
propagator measured on a = 32^^, (3 = 5.8 lattice with 
data obtained on a = 48^, (3 = 6.0 lattice (cf. Sect. V, 
Fig. 6 and Table II). 

Note that applying this procedure to several combi- 
nations of fine and coarse lattices provides us with an 
optimal scaling relation a = a(/3) for each propagator. 



APPENDIX B: TABLES 

In this appendix we present an overview of the data 
sets produced in Osaka and Berlin, the results of all 
matching fits according to Sees. V, VI, VII and Ap- 
pendix A as well as of the fits in the infrared (IR) and 
ultraviolet (UV) limits as described in Sect. IX. 



TABLE I: Lattice parameters used in this study. Configura- 
tions were generated at RCNP Osaka and HU Berlin. 
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TABLE II: Fit parameters obtained upon matching D^^ data 
from a (L",/?) = (32'', 5.8) and a (48'*,6.0) lattice for two 
different a-cuts (see Fig. 6). For comparison we also show 
the lattice spacing ratios according to Eq. (18). 



\pia\ < a 



a(5.8) 
a(6.0) 



'(5.8) 



= (6.0) 



a = 0.6 
a = 0.5 



lAltt 



1.46 
1.46 



0.989^1^ 



0.9983 



3.28 
0.92 



TABLE III: Matching the transversal gluon propagator for 
five P values (see Fig. 7): shown are the ratios of lattice spac- 
ings relative to the finest one obtained either by the matching 
procedure or according to Eq. (18); the ratios of the renormal- 
ization constants and the corresponding /dof of the fit that 
accomplishes the matching, for four choices of the a-cut. The 
lattice size is 32*. 



\pia\ < 



a = 1.0 



a = 0.8 



a = 0.6 



a = 0.5 



/3 


aifi) 
a(6.2) 






xVdof 


5.8 


i.5iil 


2.01 




5.56 


5.9 


1-31+5 


1.65 


0.9901^8 


3.83 


6.0 


1.19+^ 


1.38 


0.9931^6 


2.39 


6.1 


1 12+^" 


1.17 


0.984^3^*^ 


6.30 


5.8 


im+l 


2.01 


0.96llio 


2.16 


5.9 


1.4311 


1.65 


0.9701^5 


2.98 


6.0 


1.27t^ 


1.38 


0.973t^^ 
0.968lg 


2.61 


6.1 


1.18l«3 


1.17 


5.26 


5.8 


1.74+^ 


2.01 


0.974t^ 


1.41 


5.9 


1.46l^ 


1.65 


0.980tii 


1.70 


6.0 


1.28+? 


1.38 


0.983t22 


2.97 


6.1 


i.ist? 


1.17 


0.976127 


1.87 


5.8 


1.92i^ 


2.01 


0.963t^ 


0.585 


5.9 


1.58l^ 


1.65 


0.967tg 


2.49 


6.0 


1.39l^ 


1.38 


0.9701 1^ 


1.96 


6.1 


1.23l9^ 


1.17 


0.968l^^ 


1.59 



TABLE V: Matching the time-time gluon propagator for four 
/3 values on 32'' lattices (Osaka data): shown are the ratios 
of lattice spacings obtained by the matching relative to the 
finest one, for comparison also the ratios predicted by Necco- 
Sommer scaling, the ratios of the renormalization constants 
and the corresponding ^ jdof of the fit that accomplishes the 
matching, for two choices of the a-cut. 



\pia\ < a 



a(6.1) 



= (/3) 



(6.1) 



a = 0.6 



a = 0.5 



5.8 
5.9 
6.0 

5.8 
5.9 
6.0 



2.041^3 
1.07l^ 

2.601^1 

l-8ll?8 
1.08lg 



1.73 
1.42 
1.18 

1.73 
1.42 
1.18 



u.oz,! _g4 

0.42911^3 
^+1 



1.041 
0.3501 



0.509l^^ 



4.94 
13.9 
8.70 

1.45 
1.03 
6.19 



TABLE VI: Matching the time-time gluon propagator for 
three (3 values on 32'* lattices (Berlin data): shown are the 
ratios of lattice spacings obtained by the matching relative 
to the finest one, for comparison also the ratios predicted 
by Necco-Sommer scaling, the ratios of the renormalization 
constants and the corresponding /dof of the fit that ac- 
complishes the matching, for two choices of the a-cut. 



\pia\ 


< a 




a(/3) 
a(6.2) 


«~^(/3) 
a'^«(6.2) 


\ a(6.2) J 


xVdof 


a — 


0.6 


5.8 
6.0 


2.6ll?2 
l.llli 


2.01 
1.38 


0.292lf 
0.996l?^ 


15.7 
3.98 


a — 


0.5 


5.8 
6.0 


1.891° 
1.09l6 


2.01 
1.38 


0.487lf 
1.0211J 


11.0 
5.56 



TABLE VII: Matching the ghost propagator for four (3 val- 
ues on 32' lattices (Osaka data): shown are ratios of lattice 
spacings obtained by the matching relative to the finest one, 
for comparison the ratios predicted by Necco-Sommer scaling, 
the ratios of the renormalization constants and the /dof of 
the fit that accomplishes the matching, without (a = 2.0) and 
for two choices of the a-cut. 



TABLE IV: Fit parameters obtained upon matching the in- 
stantaneous time-time gluon propagator on a = 32'', 
P = 5.8 and a = 48'', /? = 6.0 lattice with two differ- 
ent a-cuts (see Fig. 8). For comparison we show also the 
lattice spacing ratio predicted by the Necco-Sommer scaling 
relation. 



\pia\ < a 



a(5.8) 
a(6.0) 



'(5.8) 



=(6.0) 



R^ 



a = 0.6 
a = 0.5 



1.961^^ 
1.87lfo 



1.46 
1.46 



0.476: 
0.502 



-48 
-55 
-42 
-40 



4.50 
0.503 



\pia\ 


< a 




a(f3) 
a(6.1) 


a~^(/3) 
aNS(6.i) 






5.8 




1.73 


a = 


2.0 


5.9 


I.39I3 


1.42 






6.0 


I.I8I2 


1.18 






5.8 


1.48+11 


1.73 


a = 


0.6 


5.9 


1.34+^ 


1.42 






6.0 


1.17« 


1.18 






5.8 


1.43l^? 


1.73 


a — 


0.5 


5.9 


1.35+1^ 


1.42 






6.0 


1.18+^° 


1.18 



Rz 



"(/3) 
a(6.1) 



O.955I3 

0.97ll^ 
0.985li 

0.998lii 

l-f-23 

-28 

0.987li^ 

0.9791I 

-)-l-31 
-29 



0.981] 



1.02: 



0.982: 



xVdof 



1.63 
0.840 
0.260 

0.219 
0.583 
0.420 

0.268 
0.854 
0.501 
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TABLE VIII: Matching the ghost propagator for three /3 val- 
ues on 32'* lattices (Berlin data): shown are ratios of lattice 
spacings obtained by the matching relative to the finest one, 
for comparison the ratios predicted by Necco-Sommer scaling, 
the ratios of the renormalization constants and the 'x^ /dof of 
the fit that accomplishes the matching, without (a = 2.0) and 
for two choices of the a-cut. 



\Pil\ 


< a 


P 


a((3) 
a(6.2) 


aNS(g 2) 


\ a((i.2)) 


xVdof 


a = 


2.0 


5.8 
6.0 


1.9lt} 
1.34il 


2.01 

1.38 


0.9441:1 
0.97611 


10.4 
54.6 


a = 


0.6 


5.8 
6.0 


1.2711 


2.01 
1.38 


0.959^3 
0.98712 


38.5 
121 


a = 


0.5 


5.8 
6.0 


1.40l^ 
1.18ti 


2.01 
1.38 


1.0811 
1.0211 


18.5 
19.3 



TABLE XII: Fitted UV parameters and x^/dof for the ghost 
dressing function (a = a{l3 = 6.0) = 0.1863r-o. 





[pa]min 


Cgh 


Acouifi 


7 


xVdof 


a = 0.5 


0.70 
0.75 


1.39(1) 
1.62(91) 


0.52(27) 
0.26(38) 


0.13(19) 
0.38(53) 


0.52 
0.91 


a = 0.6 


0.85 
0.9 


1.58(6) 
1.44(9) 


0.277(28) 
0.370(81) 


0.365(28) 
0.289(58) 


3.88 
2.81 


a = 2.0 


2.0 
2.5 
3.0 
3.5 
4.0 


1.50(1) 
1.51(1) 
1.56(2) 
1.55(4) 
1.55(3) 


0.313(2) 
0.300(5) 
0.274(12) 
0.273(31) 
0.270(23) 


0.313(1) 
0.319(2) 
0.329(6) 
0.330(14) 
0.332(10) 


3.96 
2.61 
1.53 
1.18 
0.87 



TABLE IX: Fitted UV parameters and /dof for the trans- 
verse gluon propagator. Data from RCNP Osaka, a = a{(3 — 
6.2) = 0.1354ro 



a = 0.5 
a = 0.6 



[pa],nin 






x'/dof 


0.5 


0.507(5) 


0.39(2) 


0.50 


0.6 


0.518(20) 


0.42(6) 


0.98 


0.5 


0.534(5) 


0.46(1) 


0.81 


0.6 


0.537(6) 


0.46(2) 


0.70 


0.7 


0.517(15) 


0.43(3) 


0.75 


0.8 


0.489(50) 


0.39(6) 


1.55 



TABLE X: Fitted UV parameters and x? /dof for the time- 
time dressing function Za with an a— cut for a = 0.5. a = 
a(/3 = 6.0) = 0.1863ro. 



data 


[pa] min 


C44a 




xVdof 


Osaka 


0.65 
0.60 
0.55 
0.50 


0.942(8) 
0.946(6) 
0.945(4) 
0.940(4) 


2.53(10) 
2.47(6) 
2.48(4) 
2.54(3) 


1.75 
1.19 
0.81 
2.18 


Berlin 


0.79 
0.75 
0.70 
0.65 
0.60 


0.968(16) 
0.962(10) 
0.962(9) 
0.948(5) 
0.940(4) 


2.09(22) 
2.19(11) 
2.19(10) 
2.38(5) 
2.49(4) 


0.23 
0.26 
0.18 
1.16 
2.22 



TABLE XI: Fitted IR parameters and X^ /dof for the time- 
time dressing function Zn with an a— cut for a = 0.5. a = 
a(/3 = 6.0) = 0.1863ro. 

data [pa]max d44a ^44 X^/dof 



TABLE XIII: Fitted IR parameters and x^ /dof for the ghost 
dressing function (a = a{f3 = 6.0) = 0.1863ro. 



\pia\ 


< a 


[pa] max 


dgha 


Kgh 


xVdof 


a = 


0.5 


0.29 
0.30 
0.34 


2.48(10) 
2.45(9) 
2.19(5) 


0.439(7) 
0.442(7) 
0.463(5) 


10.9 
8.55 
8.20 


a = 


0.6 


0.25 
0.27 
0.30 


2.57(10) 
2.54(9) 
2.50(9) 


0.435(6) 
0.437(6) 
0.440(6) 


0.62 
3.79 
4.32 


a = 


2.0 


0.24 
0.30 


2.57(9) 
2.52(9) 


0.434(6) 
0.437(6) 


5.64 
5.90 



0.20 2.07(14) 1.71(5) 0.057 

Osaka 0.25 1.52(4) 1.94(3) 19.5 

0.30 1.36(3) 2.05(2) 20.3 

0.20 4.25(32) 1.26(3) 18.2 

Berlin 0.25 2.75(8) 1.44(2) 21.2 

0.30 2.16(4) 1.57(1) 47.0 
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